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Abstract 

This  paper  presents  for  the  first  time  an  impedance-based  non-linear  model  for  lead-acid  batteries  that  is  applicable  in  all  operational  modes.  An 
overcharging  model  describes  the  accumulation  and  depletion  of  the  dissolved  Pb2+  ions.  This  physical  model  has  been  added  to  the  earlier  presented 
model  to  expand  the  model  validity.  To  properly  represent  the  charge  acceptance  during  dynamic  operation,  a  concept  of  “hardening  crystals”  has 
been  introduced  in  the  model.  Moreover,  a  detailed  gassing  and  oxygen  recombination  model  has  been  integrated.  A  realistic  simulation  of  the 
overcharging  behavior  is  now  possible.  The  mathematical  description  is  given  in  the  paper.  Simplifications  are  introduced  that  allow  for  an  efficient 
implementation  and  for  model  parameterization  in  the  time  domain.  A  comparison  between  experimental  data  and  simulation  results  demonstrates 
the  achieved  accuracy. 

The  model  enhancement  is  of  major  importance  to  analyze  charging  strategies  especially  in  partial-cycling  operation  with  limited  charging  time, 
e.g.  in  electrically  assisted  or  hybrid  cars  and  autonomous  power  supply  systems. 

©  2005  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

The  analysis  of  operation  strategies  (e.g.  in  electrically 
assisted  or  hybrid  cars  and  autonomous  power  supply  systems) 
needs  dynamical  simulation  models  for  lead-acid  batteries  as  an 
important  prerequisite.  Since  these  batteries  are  non-linear  and 
highly  dependent  on  divers  parameters  like  temperature,  state 
of  charge  (SOC)  and  short-term  history  [8,10],  the  implemen¬ 
tation  of  these  models  is  very  difficult  and  a  simple  modeling 
with  high  precision  is  consequently  nearly  impossible.  However, 
fast  computing  and  small  parameterization  effort  are  mainly 
asked. 

Electrochemical  impedance  spectroscopy  (EIS)  is  an  ade¬ 
quate  tool  for  the  development  and  the  parameterization  of  bat¬ 
tery  simulation  models  considering  these  aspects.  Small- signal 
excitation  and  the  evaluation  of  the  system’s  response  enable 
accurate  investigation  at  nearly  any  working  point.  The  comput¬ 
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ing  and  simulation  effort  of  an  impedance-based  model  can  be 
very  small  compared  with  physically  or  chemically  based  mod¬ 
els.  The  core  of  our  impedance-based  model  consists  of  elements 
that  are  related  to  physico-chemical  processes.  These  processes 
are  modeled  in  terms  of  inductors,  capacitors  and  resistors  with 
non-linear  dependency  on  temperature,  SOC  and  battery  current 
[1,2]. 

Similar  impedance-based  modeling  approaches  have  been 
presented  recently  by  several  groups  [11-15].  Barsoukov  et 
al.  [12]  and  Buller  [13]  demonstrated  the  applicability  of  this 
approach  to  several  battery  technologies  by  using  non-linear 
equivalent-electrical-circuit  models.  Andersson  et  al.  [14]  ana¬ 
lyzed  impedance  data  by  a  model  for  porous  electrodes  to  deduce 
physical  parameters,  namely  the  porosity  and  the  active  sur¬ 
face  of  the  electrodes.  For  state  of  charge  diagnosis  of  lead-acid 
batteries,  the  approach  is  employed  by  Salkind  et  al.  [15]  and 
Abolhasani  et  al.  [16],  however,  by  using  linearized  models 
that  strongly  limit  the  validity  range  of  the  diagnosis.  Formerly, 
Mauracher  and  Karden  [11]  also  proved  the  possibilities  of  this 
approach  using  a  highly  non-linear  model  whereby  the  validity 
range  was  limited  to  the  range  of  10-90%  SOC.  Additionally, 
the  model  was  limited  to  moderate  charging  currents. 
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Nomenclature 

c 

actual  ion  concentration  (molm-3) 

Cs 

equilibrium  ion  concentration  (molm-3) 

D 

diffusion  constant  (m2  s-1) 

Dcff 

effective  parameter  incl.  diffusion  constant  (s-1) 

£pb2+ 

overvoltage  at  negative  electrode  (V) 

^Pb2+pos 

overvoltage  at  positive  electrode  (V) 

/ 

radii  distribution  of  crystals  (m-4) 

F 

Faraday  constant  (96485.34)  (Cmol-1) 

A 

/ 

current  density  (Am-3) 

70,h2 

exchange  current  negative  electrode  (A) 

7o,o2 

exchange  current  positive  electrode  (A) 

fgas,neg 

gassing  current  negative  electrode  (A) 

^gas,pos 

gassing  current  positive  electrode  (A) 

Aecom 

recombination  current  (A) 

^hardening 

r  linear  scaling  factor  “hardening  crystals” 

K 

amount  of  discharged  (A  s) 

n 

valence  Pb2+  (2) 

^new 

total  number  of  not  hardened  crystals  (molm-3) 

N 

number  of  PbSC>4  crystals  (m-4) 

a 

rate  of  ions  per  time  due  to  electrochemical  pro¬ 

cess  (molm-3  s-1) 

^hardening  rate  of  hardened  crystals  (mol  m  3  s  1 ) 

^new 

rate  of  not  hardened  crystals  (molm-3  s-1) 

r 

radius  of  PbSC>4  crystal  (m) 

R 

gas  constant  (8.314472)  (Jmol-1  K-1) 

s 

rate  of  ions  per  time  due  to  chemical  process 

(molm-3  s-1) 

T 

temperature  (K) 

UT 

temperature  voltage  (V) 

fpbS04 

total  volume  of  PbSC>4  (m3  m-3) 

X0 

median  radius  (m) 

Z 

number  of  crystals  with  median  radius  (m-4) 

Greek  letters 

a  o2 

symmetric  coefficient 

7] 

overvoltage  (V) 

a 

dispersion  of  radii  distribution  (m) 

x hardening  time  constant  of  hardening  process  (s) 

L'ecom 

time  constant  of  recombination  process  (s) 

Vm 

molar  volume  PbS04  (48.22  x  10-6)  (m3  mol-1) 

In  order  to  enlarge  the  validity  range  of  this  modeling 
approach,  this  paper  presents  a  significant  improvement.  An 
overcharging  model  and  a  more  sophisticated  gassing  model 
have  been  employed  to  enable  the  analysis  of  the  charging  behav¬ 
ior  of  the  battery,  which  is  of  particular  importance  during  cyclic 
operation  with  limited  charging  time.  The  implementation  and 
parameterization  have  been  enabled  by  reference-electrode  mea¬ 
surements  inside  a  valve-regulated  lead-acid  (VRLA)  battery. 
Thereby,  investigations  on  the  dynamic  behavior  of  each  elec¬ 
trode  have  been  also  made  possible  and  will  be  additionally 
illustrated  in  this  paper. 


Without  an  overcharging  model,  the  model  output  can  devi¬ 
ate  significantly  from  the  measurement  during  charging  periods. 
This  limitation  is  inherent  to  the  impedance  method,  because 
the  battery  response  is  subject  to  a  Fourier  analysis  during  the 
impedance  measurement.  This  analysis  demands  a  constant  or 
at  most  a  time-dependent  (approximately)  linear  shift  of  the  de- 
response,  which  is  not  given  in  case  overcharging  effects  occur. 
Consequently,  these  effects  are  not  accessible  for  the  impedance 
method  and  are  thus  not  covered  by  the  impedance-based  model 
without  an  additional  overcharging  model.  Hence,  the  original 
model  significantly  underestimated  battery  voltage  when  apply¬ 
ing  charging  currents  at  high  states  of  charge  or  even  at  medium 
states  of  charge  if  a  high  charging  current  was  applied.  The 
model  validity  ends  where  the  processes  cannot  be  described 
simply  by  the  Butler- Volmer  process  and  the  accumulation  of 
acid  concentration  in  the  porous  active  material.  The  charge 
acceptance  of  the  battery  model  (without  overcharging  model) 
was  accordingly  much  too  high  during  voltage  limited  charging 
periods. 

As  an  additional  process,  which  is  identified  to  be  the  limit¬ 
ing  process  during  charging  at  high  currents  rates  or  high  states 
of  charge,  the  employed  overcharging  model  describes  the  con¬ 
sumption  of  Pb2+  ions.  The  Pb2+  ion  concentration  decreases 
during  charging  caused  by  the  generation  of  Pb  and  PbC>2  from 
Pb2+  ions.  A  liquid  phase  process  is  assumed  here.  The  dissolu¬ 
tion  process  of  the  PbSC>4  crystals  limits  Pb2+  ion  delivery  and 
finally  the  charge  acceptance  of  the  battery.  During  operation 
at  moderate  current  rates  and  low  state  of  charge,  the  dissolu¬ 
tion  rate  of  PbSC>4  crystals  is  high  and,  consequently,  it  is  not 
recognized  as  a  limiting  process. 

A  further  process  needs  to  be  recognized  which  becomes 
obvious  during  extended  experimental  tests  with  lead-acid  bat¬ 
teries.  If  a  battery  undergoes  shallow  cycles  around  a  fixed  state 
of  charge,  the  charge  acceptance  is  increasing  with  increased 
cycling  time.  This  leads  to  the  assumption  that  the  newly  formed 
sulfate  crystals  can  be  dissolved  during  charging  easier  than  sul¬ 
fate  crystals  that  have  been  formed  longer  time  ago.  This  has 
been  taken  into  account  in  the  model  by  using  a  concept  of 
“hardening  crystals”.  Hereby,  crystals  that  have  been  formed 
recently  have  a  higher  dissolution  rate  than  older  crystals. 

The  employment  of  the  detailed  gassing  model  enables  the 
description  of  the  gassing  processes  at  the  negative  and  positive 
electrode.  In  addition,  the  recombination  process  at  the  negative 
electrode  inside  a  VRLA  battery  can  be  simulated.  Thus,  the 
charge  efficiency  of  the  battery  can  be  calculated  by  regarding 
the  negative  and  positive  gassing  currents  with  respect  to  the 
terminal  current.  Moreover,  the  depolarization  of  the  negative 
electrode  can  be  simulated  which  is  forced  by  the  recombination 
process. 

The  enhanced  model  yields  high  precision  in  predicting  the 
electrical  battery  performance  even  if  high  charging  currents 
are  applied  to  the  battery.  The  parameterization  and  validation 
of  the  model  has  been  carried  out  for  a  JCI  Optima  Red  Top 
VRLA  battery  (44  Ah,  12V).  The  same  approach  can  be  also 
employed  for  flooded  lead-acid  batteries  or  prismatic  AGM  bat¬ 
teries.  Generally,  the  same  approach  should  be  also  applicable 
for  VRLA/Gel  batteries,  but  this  has  not  been  tested  so  far. 
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Fig.  1.  Equivalent  circuit  of  a  battery  half-cell  (left  hand);  equivalent  circuit  for  a  complete  lead-acid  battery  cell  (right  hand)  [2]. 


2.  The  impedance-based  modeling  approach 

For  a  basic  modeling  approach,  a  slightly  simplified  Ran¬ 
dles  equivalent-circuit  can  be  employed  ([8];  Fig.  1,  left  hand). 
The  equivalent  circuit  is  valid  for  both  battery  electrodes  with 
differing  parameters.  It  consists  of  an  internal  resistance  R[,  an 
inductance  L,  a  capacitance  C  and  a  non-linear  resistance  R.  L 
is  mainly  caused  by  the  metallic  connectors  between  the  termi¬ 
nals  and  the  electrodes  of  the  battery.  The  ohmic  resistance  R\ 
is  due  to  the  limited  conductance  of  the  contacts,  the  inter-cell 
connections,  the  grids,  the  active  masses  and  the  electrolyte.  R\ 
depends  on  SOC,  age  and  temperature  of  the  battery.  The  paral¬ 
lel  connection  of  the  capacitance  C  and  the  non-linear  resistance 
R  represents  the  double  layer  capacity  and  the  charge  transfer 
resistance. 

An  equivalent  circuit  for  a  complete  battery  cell  can  be 
obtained  by  a  series  connection  of  two  Randles  circuits.  A  few 
modifications  [1,2]  have  been  employed  for  the  impedance- 
based  lead-acid  battery  model,  namely  generalized  capacitive 
elements  (CPE:  constant  phase  element  [9]),  an  electrolyte  trans¬ 
port  model  (calculation  of  the  OCV  and  the  diffusion  overpo¬ 
tential)  and  a  parallel  current  path  for  gassing  as  the  main  side 
reaction. 

For  the  model  improvement  presented  in  this  paper,  the  equiv¬ 
alent  circuit  has  been  modified  significantly.  The  new  circuit 


considers  all  mentioned  physico-chemical  processes  including 
the  overcharging  limitations  and  side  reactions. 

3.  Overcharging  process  in  lead-acid  batteries 

3.1.  Processes  to  be  modeled 

The  overall  process  of  a  lead-acid  battery  can  be  subdivided 
into  three  different  reactions  (Fig.  2): 

1.  Electrochemical  process:  The  rate  depends  only  on  the  cur¬ 
rent  given  by  the  terminal  current.  During  discharging,  each 
pair  of  electrons  is  transferred  to  or  from  a  dissolved  Pb2+ 
ion  which  is  released  into  the  electrolyte. 

2.  Diffusion  process:  The  electrical  charge  of  each  Pb2+  ion 
must  be  counterbalanced  by  anions  like  SO42-  or  HSO4-. 

3.  Chemical  process:  The  PbSC>4  molecules  undergo  a  crystal¬ 
lization  process  and  form  PbSCU  crystals. 

During  charging,  the  process  is  reversed.  Proportional  to  the 
charging  current,  the  electrochemical  process  consumes  a  Pb2+ 
ion  flow.  The  diffusion  process  supports  these  ions  to  the  reaction 
site.  Caused  by  the  depletion  of  the  Pb2+  ion  concentration,  the 
PbSC>4  crystals  are  dissolving  [3,4]. 


current  (dissolved 
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Fig.  2.  Schematic  of  the  charge  and  discharge  reaction  on  a  microscopic  level  in  the  negative  electrode. 
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acid  density  [g/emA3] 


Fig.  3.  Solubility  of  sulfate  in  sulfuric  acid  for  different  acid  concentrations  and 
temperatures  [10]  using  stepwise  linear  interpolation. 


Assuming  a  battery  current  above  zero  (charging),  ions  are 
consumed  and  can  be  delivered  by  the  dissolution  process  of 
PbSC>4  crystals  only  (Fig.  3).  Thus,  the  ratio  of  the  actual  ion 
concentration  c  and  cs  (equilibrium  concentration,  function  of 
temperature  and  acid  concentration)  is  below  unity,  and  con¬ 
sequently,  according  to  basic  electrochemical  laws,  additional 
overvoltages  occur  (see  Eqs.  (1)  and  (2)). 


positive  electrode  :  ^pb2+p0s 


negative  electrode  : 


During  discharging,  Pb2+  ions  will  be  produced  by  the  elec¬ 
trochemical  reaction.  The  ion  concentration  c  will  increase  due 
to  the  limited  crystallization  process  of  the  PbSC>4  crystals. 
However,  due  to  the  logarithmic  behavior,  the  overvoltages  will 
change  significantly  slower  during  discharging  compared  with 
charging.  The  only  significant  impact  of  the  Pb2+  ion  concen¬ 
tration  during  discharging  can  be  observed  as  the  “coup  de 
fouet”  or  “Spannungssack”  because  the  chemical  process  is 
suppressed  at  the  beginning  of  a  discharge  process  if  no  sul¬ 
fate  crystal  seeds  are  present  already.  This  is  the  case,  if  the 
battery  has  been  very  well  charged  beforehand.  A  “coup  de 
fouet”  is  the  more  pronounced,  the  better  the  previous  charging 
has  been. 

According  to  the  illustrated  effects,  the  changing  of  the 
concentration  c  with  time  can  be  explained  by  the  following 
equation.  The  description  of  the  processes  and  the  mathematics 
follow  the  work  of  Kappus  [4]  and  its  adaptation  for  simulation 
models  by  Sauer  [7]. 


c  =  q  —  s  (3) 

The  rate  q  is  given  in  mol(m3s)_1  and  corresponds  to  the 
consumed  or  generated  Pb2+  ions  due  to  the  electrochemical 
process  per  time  (related  to  an  unified  volume  Vuni)  and  is  pro¬ 
portional  to  the  main  reaction  current,  s  is  the  number  of  ions 


that  are  dissolved  or  crystallized  from  the  PbSC>4  crystals  within 
the  chemical  process  step,  q  can  be  calculated  by  Eq.  (4)  with  I 
as  a  current  density  per  volume  Vuni,  F  is  the  Faraday  constant. 
The  Eqs.  (3)  and  (4)  are  valid  for  both  electrodes. 


/ 

q  =  — 

H  2  F 

s  can  be  calculated  according  to: 


d 
d  t 


An 


3  v 


m 


max 


r3  /(r,  t)  d r 


min 


where  r  is  the  radius  of  the  lead-sulfate  crystals, /(r)  is  the  radii 
distribution  function  of  the  lead-sulfate  crystals,  vm  is  the  molar 
volume  of  lead  sulfate,  and  rm in  and  rmax  define  the  range  of 
available  crystal  sizes.  Assuming  an  approximately  time  inde¬ 
pendent  radii  distribution  of  the  PbSC>4  crystals  (f(r,  t)  ^/(r)) 
for  small  time  steps,  the  equation  simplifies  to: 


47T 


.v  = 


3u 


m 


•rmax  Alt 

-7~(r3)f(r)dr  =  — 
at  v 


max 


r2f  f{r)  dr' 


min 


m 


min 


The  distribution/(r)  corresponds  to  the  number  of  crystals  of 
each  radius  per  volume  Vuni.  The  speed  of  radius  change  f  can 
be  calculated  by: 


vmD(c  -  c(r)) 
r 


with  D  as  a  diffusion  constant,  cs  the  equilibrium  concentration 
of  Pb2+  ions  in  sulfuric  acid  with  infinite  radius  r  and  c  the 
actual  concentration  of  the  Pb2+  ions  in  the  electrolyte  according 
to  Eq.  (3).  The  dependency  of  cs  on  the  crystal  radius  makes 
Eq.  (7)  very  difficult  to  solve.  However,  during  normal  charge 
and  discharge  operation  the  impact  of  the  dependency  of  cs  on 
the  crystal  radius  is  of  no  relevance,  therefore  it  is  not  taken 
into  account  here.1  Consequently,  Eq.  (6)  simplifies  to  (with 
c(r)  ^  cs): 

rr max 

s  =  AnD(c  —  cs)  /  r/(r)  dr  (8) 

J  6nin 


The  total  volume  of  PbSCU  Vpbso4  (related  to  Vuni)  depends 
on  the  SOC  of  each  electrode  and  is  consequently  known  from 
the  Ah  balance  taking  into  account  the  gassing  current.  K  corre¬ 
sponds  to  the  discharged  Ahs  (related  to  Vuni)  and  becomes  zero 
for  a  fully  charged  electrode  (no  sulfate  crystals  left). 


fpbSOa  = 


471  rmax  o  .  K 
—  I  r3f(r)dr  =  —vm 


min 


IF 


1  Nevertheless,  the  dependency  of  cs  from  r,  which  can  be  described  by  the 
equation  c(r)  =  cse(2orUm /RTr)  (assuming  the  validity  of  the  Gibbs-Thomson- 
relationship)  is  very  important  if  recrystallization  processes  during  open  circuit 
conditions  of  the  battery  should  be  taken  into  account.  Crystals  with  small  radius 
tend  to  dissolve  while  crystals  with  large  radius  tend  to  grow.  As  a  result,  the 
average  sulfate  crystal  radius  is  increasing.  This  processes  is  accelerated  with 
increasing  temperature  and  with  increasing  Pb2+  solubility  (and  therefore  low 
states  of  charge  as  shown  in  Fig.  3).  This  process  must  be  taken  into  account  for 
aging  models  but  this  is  out  of  the  scope  of  the  model  under  discussion  in  this 
paper. 
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Fig.  4.  Calculated  size  distribution  f(r)  of  lead-sulfate  crystals  during  discharg¬ 
ing  using  the  I\q  based  on  a  full  mathematical  approach  (published  in  [7]). 
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-  N 

discharging 

- ► 

— 

r/  m 


Fig.  5 .  Simplest  distribution  of  the  crystals  radii  with  constant  number  of  crystals 
N  and  SOC  dependent  radius  r. 

3.2.  Modeling  the  overcharging  process 


The  model  development  requires  some  simplifications  of  the 
overcharging  theory.  Namely,  the  size  distribution  f(r)  of  the 
lead-sulfate  crystals  has  to  be  approximated.2  Three  different 
modeling  approaches  have  been  tested  and  the  most  suitable 
assumption  with  an  easy  implementation  and  less  parameteri¬ 
zation  effort  has  been  chosen  finally.  A  realistic  distribution  is 
illustrated  in  Fig.  4  which  has  been  calculated  based  on  a  full 
mathematical  model  using  the  10  h  current  7io.  The  full  model 
unfortunately  takes  very  long  computation  time  and  the  imple¬ 
mentation  is  very  difficult. 


Eq.  (9)  simplifies  to: 


471 


K 


VpbS04  =  -f-Nr  =  —  Vm 
3  2  F 


(11) 


=>  r  = 


K  3 
nF  N4jt 


(12) 


Considering  Eqs.  (4),  (10)  and  (12),  Eq.  (3)  can  be  written 


as: 


c  = 


I  3/  K 


- D(c  —  cs)\  — um 3  x  42N2n2 

IF  V  nF 


(13) 


3.2.2.  Triangular  distribution 

To  reproduce  the  above  shown  realistic  distribution  (Fig.  6),  a 
triangular  distribution  with  a  medium  radius  xo  and  a  dispersion 
o  has  been  employed  for  modeling  the  overcharging  processes. 
Z  accords  to  the  number  of  crystals  with  medium  radius.  Three 
cases  have  to  be  distinguished: 


1.  The  medium  radius  xo  is  higher  than  the  dispersion  o.  The 
triangular  distribution/(r)  moves  with  decreasing  xo  and  con¬ 
stant  o  to  smaller  radii  during  charging  (see  Fig.  6,  left  hand). 
The  number  of  crystals  with  medium  radius  is  kept  fixed. 

2.  The  medium  radius  xo  is  smaller  than  the  dispersion  a  but 
greater  than  zero.  The  triangular  distribution  f(r)  moves  with 
decreasing  xo  to  smaller  radii  during  charging,  the  “left  side” 
of  the  triangle  is  cut  by  the  ordinate  and  the  number  of  crystals 
is  decreasing  (Fig.  6,  middle). 

3.  In  case  the  medium  radius  xo  reaches  zero  during  charging, 
the  condition  xo  =  0  and  a  constant  dispersion  is  assumed. 
Solely,  the  maximum  number  Zmax  decreases  further  during 
charging. 


For  case  1,  the  following  equations  can  be  found  regarding 
Eqs.  (8)  and  (9).  To  solve  Eq.  (15),  the  median  radius  xo  has  to 
be  determined  numerically  using  Eq.  (14). 


471 

kpbscu  =  —  Z 


oxl  +  Ict3x0 


(14) 


s  =  471  D(c  —  cq)Zxqo 


(15) 


3.2.1.  Equal  radius  of  all  crystals 

The  simplest  approximation  of  a  crystal-size  distribution  is 
the  assumption  of  a  fixed  number  of  A  crystals  with  equal  radius 
r  (f(r)  -  A).  The  radius  consequently  depends  on  the  total  PbSC>4 
volume  and,  thus,  depends  on  SOC  (r=/(A,  SOC)).  The  radius 
becomes  zero  for  a  fully  charged  electrode  and  increases  during 
discharging  (Fig.  5). 

Using  this  approximation,  Eq.  (5)  simplifies  to: 
d  (  4tt  f\  4tt  ? 

.v  =  —  (  - Nr  )  =  — Nr2r  =  4ixNrD(c  —  cs)  (10) 

d  t  \  3  vm  J  vm 


The  corresponding  equations  for  case  2  are  (a  numerical 
solver  is  needed  again): 


VpbS04  =  — 

3  o 


4tx  Z  /  (or  +  xo)' 


20 


~Tox'° 


(16) 


s  =  4ttD(c  —  co) 


0 


1  n  2  o 

-(xo  +  0)  -  -X0 
6  6 


(17) 


Case  3  gives  the  following  equations  at  which  no  numerical 
solvers  are  necessary: 


2  Generally  it  is  possible  to  solve  the  equation  with  numerical  methods.  This 
has  been  demonstrated  by  Sauer  in  [7].  However,  the  numerical  solution  of  the 
full  equation  is  very  difficult  and  consumes  much  computing  time. 


UpbS04  —  15^^  Zmax 
5  =  \nD(c  -  co)Zmax02 


(IB) 

09) 


958 


M.  Thele  et  al.  /  Journal  of  Power  Sources  158  (2006)  953-963 


Fig.  6.  Triangular  distribution  of  the  crystal  size,  shift  of  median  radius  xq  during  charging,  Zmax  decreases  if  .xq  =  0  and  a  remains  constant  (right  hand  figure). 


3.2.3.  Ramp  distribution 

At  higher  states  of  charge,  only  the  third  case  (Eqs.  (18)  and 
(19))  appears.  Consequently,  the  behavior  of  an  overcharging 
model  which  considers  only  this  case  does  not  differ  from  a 
model  using  the  complete  triangular  distribution.  The  imple¬ 
mentation  (Eq.  (20))  is  very  easy  and  the  parameterization  effort 
can  be  simplified  by  only  one  effective  parameter  Deff.  Conse¬ 
quently,  the  ramp  distribution  has  been  employed.  All  subse¬ 
quent  results  are  based  on  this  approach.  Eq.  (20)  is  derived 
from  Eqs.  (18)  and  (19),  no  empirical  assumptions  are  included. 

D 

s  =  10(c  -  Co)VpbS04^-  =  10(c  -  Co)CpbS04£>eff  (20) 

ac¬ 


tion  and  the  number  nnQW  corresponds  to  the  charge  Qc  of  the 
capacitor.  The  time  constant  r  =  RC  determines  the  hardening 
time  of  the  crystals.  The  variable  solubility  of  the  crystals  is 
considered  by  Eq.  (23)  using  a  slightly  modified  parameter  Deff 
(compare  Section  3.2.3).  Thus,  two  parameters  ( t  and  ^hardening) 
have  to  be  determined  by  simple  parameter  identification. 

^hardening  =  (1  T  (^hardening ^ new  ))Dq  ff  (23) 

4.  Secondary  reactions  in  VRLA  batteries 

4.1.  Processes  to  be  modeled 


3.3.  Concept  of  “ Hardening  Crystals  ” 

The  model  implementation  described  so  far  yields  an  over¬ 
charging  model  that  uniquely  depends  on  SOC.  However,  if 
a  battery  undergoes  shallow  cycling,  the  charge  acceptance  is 
increasing  with  increased  cycling  time  even  though  the  SOC  is 
approximately  stable.  The  PbS04  crystals  that  have  been  formed 
recently  appear  to  have  a  higher  dissolution  rate  than  older  crys¬ 
tals. 

A  model  modification  that  includes  the  effects  of  the  battery 
current  in  short-term  history  has  been  employed.  The  rate  of 
the  new  (not  yet  hardened)  crystals  qnQW  is  proportional  to  the 
difference  of  the  main  reaction  current  -/  and  the  rate  of  the 
hardening  crystals  ^hardening  (Eq.  (21)).  The  actual  number  of  the 
recently  formed  crystals  corresponds  to  the  integrated  value  of 
the  rate  qnQ w  (Eq.  (22)).  The  higher  nnQW  the  higher  is  the  charge 
acceptance  of  the  battery  whereat  nnQW  is  limited  to  positive 
values. 


1 

^new  —  0  ^  ^hardening 

Zr 

(21) 

^new  =  /  ^/new  tu 

(22) 

A  RC-element  has  been  employed  to  model  this  correlation 
(Fig.  7).  The  balancing  (Eq.  (21))  corresponds  to  the  nodal  equa- 


1  hardening 


Fig.  7.  Implementation  of  the  “hardening  crystals”  modeling  approach. 


Four  secondary  reactions  are  characteristic  for  a  lead-acid 
system: 

1 .  Oxygen  evolution  at  the  positive  electrode 

2.  Hydrogen  evolution  at  the  negative  electrode 

3.  Oxygen  reduction  at  the  negative  electrode 

4.  Grid  corrosion  at  the  positive  electrode 


The  first  two  reactions  can  both  be  described  by  the  Tafel 
equations. 


gas,pos 


/0,o2  exp 


nao2  q  \ 

UT  J 


(24) 


^gas,neg  —  4),H2  exP 


-nil  -  QfthPf 

UT 


(25) 


The  overvoltage  q  represents  the  difference  between  the  equi¬ 
librium  potential  of  the  respective  secondary  reaction  and  the 
present  potential  of  the  involved  electrode.  The  constants  are 
exchange  current  /o,o2  and  /o,h2,  number  of  transferred  ele¬ 
mentary  charges  n ,  symmetric  coefficients  ao2  and  au2  and  the 
temperature  T  {{RT)IF  ~  26  mV  if  T—  298  K). 

The  typical  changes  of  the  overpotential  per  current  decade 
are  —120  and  80  mV  for  hydrogen  and  oxygen  evolution  for 
flooded  lead-acid  batteries.  In  VRLA  batteries,  an  internal 
oxygen-cycle  characterizes  the  overcharging  and  float-charging 
situation.  Small  gas  channels  inside  the  electrolyte  allow  gas- 
phase  transport  of  oxygen  to  the  negative  electrode  where  it  is 
reduced  (3).  As  a  consequence  of  the  oxygen  reduction,  hydro¬ 
gen  evolution  is  largely  suppressed  and  the  overcharge  voltage  is 
mainly  determined  by  the  positive  electrode.  However,  a  certain 
evolution  of  hydrogen  at  the  negative  electrode  and  grid  corro¬ 
sion  at  the  positive  electrode  will  always  be  present,  e.g.  due  to 
the  gas  transport  limitation  inside  the  separator.  [1,5,6] 
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Fig.  8.  Impedance  measurements  of  positive  (left  hand)  and  negative  (right  hand)  electrode  at  70%  SOC  and  room  temperature  using  different  dc-currents. 


4.2.  Modeling  the  secondary  reactions  in  VRLA  batteries 

The  previous  model  topology  (see  Fig.  1,  right  hand)  was 
characterized  by  a  parallel  gassing  branch  which  combined  the 
hydrogen  and  oxygen  evolution. 

The  introduction  and  the  parameterization  of  a  more  detailed 
topology  (Fig.  10)  has  been  made  possible  by  reference- 
electrode  measurements  (inside  a  VRLA  battery)  which  allow 
separated  investigations  of  the  electrodes.  The  gas  evolution  can 
be  modeled  according  to  Eqs.  (21)  and  (22)  and  the  oxygen 
recombination  is  modeled  by  a  current  sink.  The  recombina¬ 
tion  current  depends  on  the  amount  of  oxygen  which  reaches 
the  negative  electrode.  Moreover,  the  maximum  recombination 
current  is  limited  due  to  the  limited  gas-transport  ability  inside 
the  AGM  separator.  The  effects  of  corrosion  are  neglected. 

A  phenomenological  modeling  approach  has  been  applied  to 
determine  the  recombination  current  /reCom  (see  Fig.  10).  There¬ 
fore,  a  simple  parallel  RC-element  has  been  employed,  which 
cumulates  the  generated  oxygen  in  terms  of  a  charged  capacitor 
C.  The  capacitor  can  be  discharged  by  in  due  to  the  parallel 
connected  resistor  R ,  the  time  constant  TreCom  =^recomQecom 
describes  the  time  for  the  gas  transport  through  the  separator. 
iR  is  taken  as  the  recombination  current  (/#  =  / recom)-  Finally, 
the  maximum  recombination  current  is  limited  to  /recom, max-  If 
Ir  >  /recom, max?  the  current  surplus  /  =  iR  /recom, max  corresponds 
to  an  excessive  amount  of  generated  oxygen  which  will  be  blown 
off  by  the  battery  valves.  This  mode  of  operation  can  occur  dur¬ 
ing  high  rate  charging  or  at  high  temperatures  and  should  be 
prevented  since  water  loss  is  caused  which  leads  to  a  reduced 
service  life  of  the  battery,  /recom, max  can  change  during  the  life¬ 
time  of  a  battery.  It  is  very  low  for  a  brand  new  battery  because 
only  a  small  number  of  gas  pores  is  available.  Within  increasing 
loss  of  water,  the  number  of  gas  channels  increases  and  therefore 
hecom,max  increases  as  well. 

5.  The  improved  battery  model 

5.1.  Model  structure 

Reference-electrode  measurements  allowed  also  separated 
impedance  measurements  on  the  electrodes  whose  results  are 


depicted  in  Fig.  8.  Both,  the  positive  electrode  and  the  negative 
electrode  are  characterized  by  a  non-linear  current-dependent 
dynamic  behavior. 

The  positive  impedance  shows  a  small  current-independent 
semi-circle  for  high  measurement  frequencies  /ac  and,  for 
increased  /ac,  a  second  significant  larger  current-dependent 
semi-circle.  The  measured  non-linear  behavior  of  the  positive 
electrode  is  illustrated  in  Fig.  9  (left  hand,  “Rp0s_2  (spectra)”) 
and  can  be  explained  by  a  Butler- Volmer  behavior.  Also  given 
is  the  current- dependent  behavior  of  the  first  semi-circle  of  the 
negative  electrode.  It  corresponds  to  an  already  observed  char¬ 
acteristic  [1]  and  can  be  modeled  by  a  phenomenologically 
assumed  tanh-function.  The  behavior  of  the  negative  electrode 
at  lower/ac  is  determined  by  a  second  “circle”  which  reaches  the 
positive  imaginary  region.  This  part  of  the  spectrum  has  been 
ignored  in  the  parameterization  of  the  model. 

According  to  the  spectra  and  the  above  mentioned  pro¬ 
cesses,  a  complex  model  structure  has  been  developed  (Fig.  10). 
The  electrodes  are  completely  separated  from  each  other.  The 
dynamic  behavior  of  both  electrodes  is  considered  by  RC  ele¬ 
ments  with  partly  non-linear  resistors.  Actually,  the  quiescent 
electrode  potential  should  be  connected  in  parallel  to  the  dou¬ 
ble  layer  capacitance  and  in  series  with  the  charge  transfer 
resistance.  Since  a  clear  identification  of  both  effects  is  not  pos¬ 
sible  (each  electrode  provides  two  RC  elements),  the  quiescent 
potential  has  been  connected  in  series  with  the  RC  elements. 
Nevertheless,  the  dynamic  behavior  is  not  influenced  by  this  sim¬ 
plification.  Also,  the  overcharging  models  have  been  employed 
by  two  voltage  sources  £Pb2+pos  and  £pb2+pos  whose  values  are 
calculated  according  to  Eqs.  (1)  and  (2).  The  secondary  reactions 
are  considered  as  illustrated  above. 

5.2.  Parameterization  of  the  model 

Detailed  information  on  the  impedance-based  parameteriza¬ 
tion  procedure  of  the  RC  elements  and  the  time-domain  param¬ 
eterization  of  the  electrolyte  transport  model  can  be  found  in 
the  literature  [1,2].  The  parameterization  of  the  RC  elements  is 
based  on  galvanostatic  impedance  measurements  using  several 
bias  dc-currents  at  different  SOC  and  temperatures.  Thus,  the 
non-linear  dynamic  behavior  of  the  battery  can  be  determined. 
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current  /  A 


Fig.  9.  Non-linear  behavior  of  the  electrodes  (70%  SOC  at  room  temperature);  “spectra”-values  are  extracted  from  measured  impedance  spectra,  “modeF’-curves 
correspond  to  a  Butler- Volmer  behavior  (positive  electrode)  and  a  phenomenologically  modeled  behavior  (negative  electrode). 
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Fig.  10.  The  complete  battery  model  topology  with  separated  electrodes. 


The  parameterization  of  the  transport  model  consists  of  the  eval¬ 
uation  of  defined  discharging  steps  at  several  SOC. 

An  effective  parameter  Deff  could  be  derived  for  a  ramp  dis¬ 
tribution  of  the  sulfate  crystals  size  (Eq.  (20)).  This  parameter 
has  been  parameterized  for  different  states  of  charge  (60,  70  and 
90%)  using  alternating  charging  and  discharging  steps.  The  cur¬ 


rent  value  has  been  increased  until  overcharging  effects  could 
be  identified,  respectively,  the  maximum  accepted  charging  cur¬ 
rent  has  been  found  (see  Fig.  11,  left  hand).  The  parameter  Deff 
is  given  in  Fig.  1 1  (right  hand).  Since  only  the  values  at  60,  70 
and  90%  SOC  have  been  determined,  the  remaining  values  are 
extrapolated. 


Fig.  11.  Cell  voltage  during  alternating  charging  and  discharging  steps  at  90%  SOC  with  10  s  duration,  begin  of  overcharging  effects  at  8.8  A  (4  x  /20)  (left  hand); 
effective  parameter  De g-  determined  at  60,  70  and  90%  SOC,  all  remaining  values  are  linearly  inter-  and  extrapolated  (right  hand). 
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For  the  parameterization  of  the  gassing  model,  different  con¬ 
stant  currents  have  been  applied  to  the  fully  charged  battery.  The 
resulting  potentials  of  the  single  electrodes  have  been  measured 
after  a  minimum  of  a  24  h  period  whilst  applying  a  constant 
current  (Fig.  12).  The  positive  electrode  shows  an  almost  expo¬ 
nential  behavior,  whereas  the  negative  electrode  is  characterized 
by  two  different  effects  dependent  on  the  level  of  the  applied  cur¬ 
rent.  For  lower  currents,  sufficient  oxygen  can  be  transported 
to  the  negative  electrode  and  the  applied  current  is  consumed 
for  oxygen  reduction.  The  electrode  is  depolarized.  For  higher 
currents,  the  applied  current  is  not  only  consumed  for  oxygen 
reduction  but  also  for  hydrogen  generation.  The  negative  elec¬ 
trode  polarizes.  The  parameterized  modeling  of  the  secondary 
reactions  is  given  with  solid  lines.  The  generally  known  slopes 
of  oxygen  and  hydrogen  generation  are  depicted  with  dashed 
lines.  Further  information  can  be  found  in  [17]. 

6.  Simulation  results 


-1.0  -0.8  -0.6  -0.4  -0.2  0.0  1.6  1.7  1.8  1.9  2.0 

overpotential  /  V 


Fig.  12.  Overpotential  of  negative  and  positive  electrode  related  to  hydrogen 
reference-electrode;  solid  lines:  modeling  of  gassing  behavior  (VRLA  battery), 
dashed  lines:  Tafel  slopes  according  to  literature  [10]. 


For  verification,  a  dynamic  current  profile  with  alternat¬ 
ing  pulses  (10  s  duration)  and  rest  periods  in  between  has 
been  applied.  Values  from  ±2.2  to  ±61.6  A  have  been  used 
which  correspond  to  multiples  of  the  nominal  battery  current 


/n  =  ho  =  2.2  A  (1  x  /20  to  28  x  /20).  The  verification  has  been 
carried  out  using  a  voltage  limitation  at  2.43  V/cell.  A  PID  con¬ 
troller  has  been  added  to  the  simulation  model,  so  that  the  battery 
current  has  been  controlled  during  the  voltage  limited  charging 


Fig.  13.  Simulation  and  measurement  data  at  90%  SOC,  repetition  of  several  current  steps  with  ±2.2  A  up  to  ±61.6  A  current  value,  charging  periods  are  voltage 
limited  at  2.43  V  (red-top  Optima  battery,  12  V,  44  Ah). 


Fig.  14.  Left  hand:  measurement  and  simulation  of  battery  current  during  constant  voltage  charging  without  “hardening  crystals”  concept  corresponding  to  Fig.  13; 
right  hand:  comparison  of  SOC  calculations  (“measurement”:  integration  of  the  measured  terminal  current  (red-top  Optima  battery,  12  V,  44  Ah)). 
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Fig.  15.  Measurement  and  simulation  results  at  70%  SOC  and  room  temperature  (red-top  Optima  battery,  12  V,  44  Ah),  repetition  of  current  steps  with  ±30.8  and 
±61.6  A  current  value,  charging  periods  are  voltage  limited  at  2.43  V  (left  hand);  battery  current  during  the  constant  voltage  charging  periods  without  “hardening 
crystals”  concept  (right  hand). 


Fig.  16.  Left  hand:  simulated  and  measured  current  during  constant  voltage  charging  corresponding  to  Fig.  14  but  using  the  “hardening  crystals”  modeling  approach; 
right  hand:  significantly  improved  SOC  calculation. 


periods.  This  procedure  is  very  similar  to  battery  charging  equip¬ 
ment  during  constant  voltage  charging.  Thus,  the  simulated  and 
measured  voltages  during  the  limitation  are  equal  and,  for  veri¬ 
fication  of  the  model,  the  controlled  currents  have  consequently 
also  to  be  compared  during  the  voltage  limitation. 

Initially,  the  measured  and  the  simulated  voltage  responses 
are  compared  at  90%  SOC  (Fig.  13).  In  spite  of  the  very  high 
resolution,  only  slight  deviations  become  obvious  (<10mV). 
The  left  hand  figure  illustrates  the  first  two  different  current 
values  (±2.2  and  ±4.4  A)  whereat  the  voltage  limitation  is  not 
reached.  The  figure  on  the  right  hand  gives  the  results  of  the 
last  three  current  values  (±15.4,  ±30.8  and  ±61.6  A).  Both,  the 
measurement  and  the  simulation  reach  the  voltage  limitation. 

Moreover,  the  controlled  currents  are  compared  in  Fig.  14 
(left  hand).  Without  modeling  the  mentioned  “hardening  crys¬ 
tals”  effect,  the  overcharging  model  is  only  dependent  on  SOC. 
Thus,  the  simulated  current  shows  no  dependency  on  the  bat¬ 
tery  current  history  (e.g.  no  dynamic  behavior  can  be  observed 
during  the  10  s  charging  steps).  Nevertheless,  the  use  of  an  over¬ 
charging  model  (even  though  dependent  on  SOC  only)  limits 
the  charge  acceptance  at  all  and  improves  the  SOC  calculation 
consequently  (right  hand). 


The  same  verification  profile  has  been  applied  to  the  battery 
at  70%  SOC.  Again,  only  negligible  deviations  between  the  sim¬ 
ulated  and  the  measured  voltages  are  observable  (Fig.  15).  The 
deviations  of  the  battery  currents  during  the  constant  voltage 
charging  periods  (again  without  using  the  concept  of  “harden¬ 
ing  crystals”)  are  smaller,  but  qualitatively  the  same  as  detected 
at  90%  SOC. 

Finally,  the  current  deviations  can  be  minimized  by  the 
employment  of  the  “hardening  crystals”  approach.  Fig.  16  illus¬ 
trates  the  results  at  90%  SOC  corresponding  to  Fig.  14.  Both,  the 
absolute  value  of  the  accepted  charging  current,  which  grows 
with  increased  discharging  current  in  short-term  history,  and 
the  dynamic  behavior  during  the  charging  steps  can  be  simu¬ 
lated  nearly  analog  to  measurement  data.  The  SOC  calculation 
is  improved  accordingly. 

7.  Conclusions  and  summary 

A  significantly  improved  modeling  approach  for  lead-acid 
batteries  has  been  illustrated.  Primarily,  the  overcharging  pro¬ 
cesses  have  been  discussed  and  useful  mathematical  approaches 
have  been  given  for  modeling.  Finally,  the  concept  of  “harden- 
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in g  crystals”  has  been  identified  as  an  important  prerequisite  for 
the  simulation  of  the  batteries  charge  acceptance.  The  concept 
validation  has  been  performed  by  using  measurements  on  a 
VRLA/AGM  battery. 

However,  more  theoretical  explanation  and  experimental  evi¬ 
dence  for  the  “hardening  crystal”  concept  would  be  beneficial 
for  the  understanding  and  adequate  parameterization.  Moreover, 
simplifications  regarding  the  model  topology  and  the  parameter¬ 
ization  process  are  imaginable  particularly  for  flooded  lead-acid 
batteries.  Since  the  recombination  process  is  negligible  there, 
the  separation  of  both  electrodes  is  not  absolutely  necessary  any¬ 
more.  The  gassing  and  a  overcharging  model  for  both  electrodes 
could  be  merged,  removing  the  need  for  reference-electrode 
measurements.  The  simplified  model  structure  (Fig.  1,  right 
hand)  could  be  reused. 
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